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INTRODUCTION 

Several topics relating to combustion limits in premixed flames at reduced gravity have been 
studied. These topics include: (1) flame balls; (2) numerical simulation of flame ball and planar 
flame structure and stability; (3) experimental simulation of buoyancy effects in premixed flames 
using aqueous autocatalytic reactions; and (4) premixed flame propagation in Hele-Shaw cells. 

STRUCTURE OF FLAME BALLS AT LOW LEWIS-NUMBER (SOFBALL) 

Successful flame ball experiments were conducted on the STS-83 and STS-94 Space Shuttle 
missions in 1997. Among the unexpected results were (1) flame balls survived much longer than 
expected based on pre-flight estimates - in most cases the entire 500 second test duration; (2) when 
multiple flame balls were present, they drifted apart from each other (although their “center of 
mass” was nearly motionless in many cases); (3) flame balls were very sensitive to small 
accelerations (typically 100 pg for 1 second) resulting from Orbiter vernier reaction control system 
(VRCS) thruster firings; and (4) remarkably, all flame balls in all mixtures tested produced 
between 1 and 2 Watts per ball of radiant power per ball. 

The drift of adjacent flame balls was attributed to the enthalpy gradient each ball imposes on its 
neighbors. A model was developed by the PI in conjunction with J. D. Buckmaster. The slope of 
separation vs. time plots (e.g. Fig. 1), are close to the theoretically prediction of 1/3, and the 
quantitative agreement is reasonable considering the simplicity of the model. The mean and 
standard deviation of the slope for all relevant tests is 0.334 with a standard deviation of 0.060. 

The flame balls were found to respond ballistically to VRCS impulses (Fig. 2). The impulse 
(change in velocity) imparted to the ball is about the same as the acceleration impulse (integral of 
acceleration over time). The ball velocity change is somewhat greater (by a factor of about 2) than 
the acceleration impulse. The maximum possible factor for a spherical bubble of fluid having very 
low density compared to the surrounding fluid is 2 because of the “added mass” effect. This 
change in velocity then decays on a time scale of tens of seconds, which is comparable to the 
viscous time scale associated with the flame ball and its surrounding hot gas field. 



Time (seconds) 


Figure 1 . Examples of flame ball interactions resulting Figure 2. Effect of impulsive orbiter accelerations on 
in mutual repulsion. flame ball drift velocity. 

NUMERICAL SIMULATION OF FLAME BALL AND PLANAR FLAMES 

In order to predict the properties of flame balls, a one dimensional, time-dependent flame 
code with detailed chemical and transport sub-models was employed to solve the usual 
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nonsteady equations for energy and species conservation in spherical geometry at constant 
pressure. It was found that the predicted flame ball radii using different chemical models vary 
widely, even though all of these models predict the burning velocities of steady plane flames 
quite accurately. Also, reabsorption of emitted radiation was found to be a dominant effect in 
mixtures diluted with CO2 or SF 6 , in that calculations assuming infinite absorption coefficient for 
the diluent showed much better agreement with experiment that calculations assuming optically- 
thin diluent radiation (Fig. 3 ). 

A numerical study of premixed-gas flames in mixtures of CH4, O2, N2 and CO2 using detailed 
chemical and radiative emission-absorption models was conducted in conjunction with Prof. 
Yiguang Ju of Tohoku University in Sendai, Japan. It was found that reabsorption of emitted 
radiation led to substantially higher burning velocities and wider extinction limits than 
calculations using optically-thin radiation models, particularly when CO2, a strong absorber, is 
present in the unbumed gas (Fig. 4 ). Two heat loss mechanisms that lead to flammability limits 
even with reabsorption were identified: (1) differences in the absorption spectra of the reactant 
and product molecules and (2) the broadening of emission spectra at flame temperatures 
compared to ambient temperature. Via both mechanisms some net upstream heat loss due to 
radiation will always occur, leading to extinction of sufficiently weak mixtures. It is concluded 
that fundamental flammability limits can exist due to radiative heat loss, but these limits are 
strongly dependent on the emission-absorption spectra of the reactant and product gases and 
their temperature dependence, and cannot be predicted using gray-gas or optically-thin models. 



Figure 3. Predicted flame ball radii in H 2 -O 2 -CO 2 
mixtures (H 2 :02 = 1:2) including and excluding CO 2 
radiation, along with measured flame ball radii from 
aircraft fig experiments, JAMIC drop-tower 
experiments, and the STS -83/94 space experiments. 


Figure 4. Predicted values of burning velocity and peak 
flame temperature in CH 4 - (0.21 0 2 + 0.49 N 2 + 0.30 
CO 2 ) mixtures under adiabatic conditions, with 
optically-thin radiative losses, and including 
reabsorption effects. 


“LIQUID FLAMES” 

With NASA Code UG support, the PI has introduced the use of aqueous autocatalytic 
chemical reaction fronts for the experimental simulation of combustion processes. Theses fronts 
exhibit very little density change across the front, have simple chemistry, are unaffected by heat 
losses (since the front is nearly isothermal) and have high Schmidt numbers, allowing the front to 
remain “flamelet-like” even in the presence of very strong flow disturbances or turbulence. Thus, 
such fronts are very useful for experimental study of combustion under conditions far more 
readily simulated by available theoretical and numerical models. 

Hele-Shaw cells are frequently employed to study buoyancy effects in fluid systems in a 
simple quasi two-dimensional geometry where the flow is governed by a linear equation (Darcy’s 
law). Thus, autocatalytic reactions in Hele-Shaw cells represent the simplest possible 






experimental realization of the interaction of a propagating front with buoyancy-induced 
convection. Fingering-type instabilities were observed (Fig. 5), which is very surprising since the 
classical Saffman-Taylor mechanism is not present (since there is no viscosity difference across 
the front). The fingering wavelength is independent of the front propagation rate (Sl) and the cell 
thickness (w) and is inversely proportional to the square root of the cosine of the angle of the cell 
from vertical. The only viable explanation of this behavior is a surface tension at the interface 
whose magnitude is about 0.005 dyne/cm - about 14.000 times smaller than a water-air interface. 

An estimate of the magnitude of the buoyancy-induced flow disturbances was derived based 
on the Saffman-Taylor model, resulting in u’/Sl *» (7t/12)g8w 2 /vSL where u’ is the effective 
“turbulence intensity”, g the gravitational acceleration, 8 the fractional density change, and v the 
kinematic viscosity. The “turbulent burning velocity” (St) is plotted as a function of U in Fig. 6. 
These results show that the Yakhot model Sj/Sl = exp ((u7Sl)/(S-j/Sl)) 2 fit the experimental data 
for this experiment, as it did our previous results from several different forced-turbulence flows. 
These results suggest a rather simple description of the role of buoyancy on the front 
propagation in this simple chemical and hydrodynamic system. 



Figure 5. Upward propagating autocatalytic front in a Figure 6. Correlation of wrinkled front speed (Ut) with 
Hele-Shaw cell. Cell thickness (w) =1.0 mm, Sl = buoyancy parameter (Ub). 

0.17 mm/s. Upper image: 10 seconds after initiation; 
lower image: after reaching quasi-steady propagation 
condition. Width of cell is 200 mm. 

PREMIXED-GAS FLAME PROPAGATION IN HELE-SHAW CELLS 

As a complement to the experiments on chemical fronts in Hele-Shaw cells, premixed-gas 
flames in Hele-Shaw cells were also examined. Significantly, wrinkling was observed even for 
downward propagating (buoyantly stable) flames and flames having high Lewis number 
(diffusive-thermally stable) (Fig. 7a). The burning rates (S T ) of these flames are quite different 
from their laminar, unwrinkled values (Sl). Values of St/Sl in the quasi-steady stage were higher 
for upward vs. downward propagation, but only weakly dependent on Lewis and Peclet numbers 
(Fig. 7b). The first of these effects is consistent with the Joulin-Sivashinsky (JS) model, whereas 
the last is contrary to the JS predictions concerning the effects of heat losses. These results 
show that even for diffusively stable mixtures, at microgravity thermal expansion and viscosity 
changes across the front will lead to flame instabilities. These results also indicate that the 
behavior of flame propagation in narrow channels such as crevice volumes in premixed-charge 
internal combustion engines (the source of most unbumed hydrocarbon emissions) may be quite 
different from that inferred from simple laminar flame experiments. 
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(a) (b) 

Figure 7. Characteristics of flames in Hele-Shaw cells, (a) direct images (flames propagate from left to right) Cell 
width (vertical direction in these images) 39 cm. Cell length (horizontal direction in these images) 60 cm, but 
images are cropped to show only flame front Images from left to right 7.2% CH 4 in air, horizontal propagation; 
7.1% CH< in air, upward propagation; 7.1% CH 4 in air, downward propagation; 3.0% CjHs in air, horizontal 
propagation, (b) Correlation of wrinkled front speed (St/Sl) with Peelet number. 
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